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^ ■ Abstract 
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Neutrino oscillations with parameters Attiq = (2 — 20) • 10 ^ eV"^, sin"^ 20q > 0.65, 

m ; 

CN i relevant for large mixing MSW solution of the solar neutrino problem can lead to an 

O '. 

On ' observable (up to 10 - 12 %) excess of the e-like events in the sub-GeV atmospheric 



ON 



neutrino sample. The excess has a weak zenith angle dependence in the low energy 

part of the sample and strong zenith angle dependence in the high energy part. The 
D , excess rapidly decreases with energy of neutrinos, it is suppressed in the multi-GeV 

sample. These signatures allow one to disentangle the effect of the oscillations due to 
K> , solar Am? from other possible explanations of the excess. The up-down asymmetry of 

C^ ' the excess may change the sign with energy being positive in the sub-GeV region and 

negative in the multi-GeV range. Predicted properties of the excess are in agreement 

with the SuperKamiokande data. 



1 Introduction 

The SuperKamiokande (SK) collaboration [|I|, ^ continues to accumulate the data which 
strengthen the evidence of the muon neutrino oscillations. Whole variety of the data can 
be fitted well assuming z/^ ^^ v^ channel with the maximal or close to maximal mixing: 

Am^,„ = (1.5 - 8) X 10^^^ , sin^ 2^,,^ > 0.8. (1) 

At the same time, there are some facts which indicate that oscillations of the atmospheric 
neutrino is not reduced completely to two neutrino case and the electron neutrino is involved 
in the oscillations too. The v^ ^->- z/g oscillations as the sub-leading mode are possible and 
probably desired. 

• The data shows an excess of the e-like events in the sub-GeV sample. The SK col- 
laboration accounts the excess by up-scaling the overall normalization of the neutrino 
fluxes. However, recent cosmic ray measurements [§] indicate that possible increase of 
the neutrino fluxes was probably overestimated. Moreover, it seems that the excess 
depends on the neutrino energy (which is impossible to explain by overall normaliza- 
tion). The largest excess is in the low energy part of the sample {p < 0.4 GeV) and it 
is smaller in the high energy part. 

• There are some indications of the up-down asymmetry of the e-like events which 
changes with energy. The asymmetry is positive in the low momenta region and it 
tends to be negative at high energies (multi-GeV range). 

• The data show strong zenith angle dependence of the /i-like events which imply rel- 
atively small oscillation effect for down-going neutrinos. In the two neutrino frame- 
work it is difficult to reconcile this fact with rather low value of the double ratio 
i?^/e = (yu/e)/(/i/e)A/c integrated over zenith angle [^]. (Although recent data give 
higher -R^/e the problem is not completely settled down.) In the case of three neutrino 
oscillations the i/^-signal can be further suppressed and the i/g-signal enhanced. 



There have been a number of studies of the atmospheric neutrino oscillations with the 
three (or more) neutrino mixing [^, |^, |^, ||, ^, |1^, |lT|. Most of them where performed 
in the framework of the so called "one level dominating scheme" when the mass splitting 
between two lightest states, Am2i, is neglected. In this case the oscillations in the sub- 
leading channel occur due to mixing of the electron neutrino in the heaviest mass state 
with mass splitting A?7i^^^. The effects of the sub-leading channel are restricted by the 
CHOOZ result |12|. They are reduced to vacuum oscillations for the sub-GeV sample. For 
the multi-GeV sample the Earth matter effect becomes important which may enhance or 
suppress the oscillations PH]. 



In several papers all mass splittings have been taken into account. In Ref. [|TT| generic 
3z/ effects have been considered with relatively large Amli for lightest states (outside the 
region of solar neutrino solutions). In Ref. [l^ the attempts have been made to explain 
the zenith angle dependence of the /i-like events by the z/^ ^->- z/g oscillations with small 
Am^i ~ 10~^ eV^ and large mixing. The channel z/^ ^^ u^ with large ArTigg ~ 0.2 eV^ 
implied by the LSND result leads to overall suppression of the signal. It was marked that 
for Am2i < 10~^ eV^ the matter of the Earth significantly suppresses the oscillation effect. 
In this scheme one expects also strong zenith angle dependence of the e-like events which 
contradicts the SK data. 



It was marked in Ref. WM that the effect of the sub-leading oscillations driven by AttIq 
responsible for the solar neutrino deficit is significant only for the sub-GeV events and the 
size of effect at the level of the statistical errors. In Ref. |]15[ it was argued that the excess 



of the e-like events in the sub-GeV sample favors large mixing MSW solution of the solar 



neutrino problem. However analysis of 535 days of the SuperKamiokande data in Ref. [|16 
leads to conclusion that there is no difference between the large and small mixing solutions. 

The impact of the solar neutrino sector on the atmospheric neutrino fluxes has been 
also studied in the context of exact parity model H^ . 



In this paper we will study in details the effects of oscillations of the atmospheric neutri- 
nos with parameters (Amg, sin^ 2^©) of the large mixing MSW solution of the solar neutrino 
problem. The analysis of recent data (including the spectral and zenith angle information) 
leads to the following values 0, |18 



Am^ = (2 - 20) ■ 10"^ eV^ , sin^ 26^ = 0.65 - 0.95 . (2) 

This region corresponds to two neutrino mixing. It is slightly modified if there is an admix- 
ture of z/g flavor in the third state which satisfies the CHOOZ bound ||12|: |f/e3p = sin^ 6*13 < 
0.05. 

There are two motivations of present study: 

(1). Recent SK data give some indications in favor of the large mixing MSW solu- 
tion (fiat distortion of the recoil electron spectrum, day-night asymmetry, fiat zenith angle 
distribution of the night signal, etc.) [0, |19 . 



(2). Studies of the atmospheric neutrinos enter now new stage of precision measurements 
when the data become sensitive to sub-leading effects. 

The paper is organized as follows. In Sec. 2 we find general expressions for the atmo- 
spheric neutrino fluxes in presence of the three neutrino oscillations. In Sec. 3 we discuss 
properties of relevant two neutrino probabilities. In sect. 4 we calculate the excess of the 
e-like events for different energy ranges and study its properties. In sect. 5 we consider 
influence of the sub-leading oscillations on the //—like events and -R^/e- In Sec. 6 we discuss 
the results and perspectives. 

2 Three-flavor oscillations of atmospheric neutrinos 

We consider the three-flavor neutrino system with hierarchical mass squared differences: 
Am^i = AuIq « Am3]^ = Am^^^ (see Eqs. (|l],|D). The evolution of the neutrino vector of 



state Uf 



is described by the equation 



.duf _ (UM'^W 



+ V\vf, 



(3) 



dt \ 2E 

where E is the neutrino energy and M^ = diag{0, Amli, Am^^) is the diagonal matrix 
of neutrino mass squared eigenvalues. V = diag{Ve, 0, 0) is the matrix of matter-induced 
neutrino potentials with Ve = ylGpN^^ Gp and A'"e being the Fermi constant and the 
electron number density, respectively. The mixing matrix U is defined through vj = Uum, 
where Um = {^i, ^2, ^3)^ is the vector of neutrino mass eigenstates. It can be parameterized 
as f/ = U2zUi'^Ui2. The matrix U-ij = Uij{6ij) performs the rotation in the ij- plane by the 
angle 6ij. We have neglected possible CP- violation effects in the lepton sector which are 
suppressed in the case of the mass hierarchy. 

Let us introduce new states z/ = (z/g, z/2, t's)"^ obtained by performing the f/23 - rotation: 
Uf = L^23^'- The Hamiltonian H that describes the evolution of the vector u can be found 
from Eq. (|^): 

H = l^U..U,oM^Ulul + V . 
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Let us assume that mixing of the electron neutrino in the heavy state is negligiblej^ so 
that f/i3 ^ 1. In this case we get explicitly 
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(ci2 = cos6'i2, S12 = sin6'i2, etc.) According to Eq. (^, the z/3 state decouples from the rest 
of the system and evolves independently. Therefore the S-matrix (the matrix of amplitudes) 
in the basis (t'e, ^^2, ^^3) has the following form : 
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(5) 



^ We comment of the effect of this mixing in Section 6. Results of the detailed studies will be published 
elsewhere pOl. 



where 



A. 



33 = exp(-«03J 
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and L is the total distance traveled by the neutrinos. The {Uf., V2) subsystem evolves ac- 
cording to the 2x2 Hamiltonian {v^. — v^ sub-matrix in Eq. (^)). The latter depends on the 
potential V^, mixing angle Qyi and the mass squared difference Am^i. Let us denote by 
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the probability of the z/g <-> z/2 oscillations. For antineutrinos we have P2 = -P2(^K)- 

The S'-matrix in the flavor basis can be obtained from Eq. (^) by t/23-rotation: U23SU23. 
It gives the probabilities of flavor oscillations as Pi^a -^ ^p) = |(f^23'S'?723)a/3p- The appro- 
priate probabilities equal 
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where 



is the phase difference of the amplitudes ^22 and ^433: 02 = argA22 is the phase of ^22; 03 
is defined in Eq. (|^). 

Using the probabilities given in Eqs. (18|-|T0D one can find modifications of the atmospheric 



neutrino fluxes due to the oscillations. Let F^ and -F° be the electron and muon neutrino 
fluxes at the detector in the absence of oscillations. Then the fluxes in the presence of 
oscillations can be written as 



I + P2 (re 3-1 
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2 2^23(l-\/l-Acos( 



'11^ 



(12) 



where 



is the ratio of the original muon and electron neutrino fluxes. Here Q^ is the neutrino zenith 
angle. 



For antineutrinos P2 should be substituted by P2 in Eqs (|| - |I2]). The solution of the 
solar neutrino problem implies that the resonance is in the neutrino channel, therefore the 
mixing and the transition probability are smaller for antineutrinos: P2 > P2- 

As follows from Eq. (|11]) the effect of oscillations on the electron neutrino flux is pro- 
portional to the factor {rc^^ — 1). Therefore one can have either an excess or a deficiency of 
the e-like events depending on values of r and C23. The ratio r depends both on the zenith 
angle and on the neutrino energy. For r = 2 which corresponds to the sub-GeV sample, 
there will be an excess of e-like events for 6*23 < 45° and a deficiency for ^23 > 45°. The SK 
best fit is ^23 = 45°; in this case there would be no deviation from the prediction for r = 2. 
In the multi-GeV range r is typically 3 - 3.5 in the vertical direction with averaged over 
zenith angle value about 2.5. 

3 The transition probability 

We have calculated the two neutrino transition probabilities P2, P2 using the parameter- 
ization of the distribution of density in the Earth from Ref. [^. For the analysis of the 
results, it is worthwhile to consider transitions of neutrinos in the Earth as oscillations in 
medium which consists of several layers with constant densities. This, in fact, gives rather 
good approximation to the exact results. 

The depth and the length of oscillations are determined by 

sin^ 2^^ = sin^ 2^12 (^^^ j , ^™ = ^^^ (13) 



where 

AH 



Arriai , . \ / „ Amh ^ 



\\y^os2e,2^-Kj +I^sm29,,^j (14) 

is the level splitting (difference between the eigenvalues of H). The resonance energy equals 

Am^iCos26'i2 ^ / Am^ \ /'2.0q/cm^\ ^ , , 

Er = ^^— = 0.17 GeV 1^ — ^^ cos 29u. 15 

In the mantle, for typical value A?n2i = 5 ■ 10~^ eV^ and for sin^ 29i2 = 0.8 we get 
Eji = 0.08 GeV which is below the threshold of the sub-GeV range. Therefore for 
Am\^ < 5 ■ 10~^ eV^ the oscillations are in the matter dominated regime when the po- 
tential is substantially larger than the kinetic term: V ^ Aw? J2E. In this case, the depth 
of oscillations is roughly proportional to (Arra^)^, and the oscillation length, Z^, is close to 
the refraction length, /q? and only weakly depends on the energy: 



sin^ 2^„ ~ sin^ 2^12 I — i^ 1 , /„ ~ /q = ^. (16) 
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For the multi-GeV range the approximation of Eq. (|T6D works for Am^ as big as 10 ^ eV^. 

We find that the maximal neutrino oscillation effect in the mantle is achieved at cos Qy ~ 
—0.35 and the effect is zero at cosO,y ~ —0.6. For cosO,y < —0.84 neutrinos cross both the 
mantle and the core of the Earth. The interplay of the oscillations in the mantle and in the 
core leads to some enhancement of the transition probability in spite of larger density of the 
core. The oscillation effects in the antineutrino channel are smaller by factor 2-3. 

The expressions of the Eq. (|16]) are valid for small sin^ 2Qm- With increase of Am} ^ the 
increase of sin^ 2^^, and consequently, the probabilities is slowing down. In the neutrino 
channel the depth approaches one in the resonance. In the antineutrino channel the mixing 
also increases but it is always below vacuum mixing. 

4 Excess of the e-like events 

In what follows we will calculate the dependences of the excess of e-like events on the zenith 
angle of electron, Og. The general expression for the number of e-like events, A^e as a function 



of 6p is 



N,ocJ2j dE,dEJ{cose,)dh FeiE,,Q,)^^ (17) 

vv 

^{Qe,Qu,E,)Ke{h,cose„E,)e{E,) , (18) 

where Eg is the atmospheric z/e-flux at the detector given in Eq. ([Tl|) (the fluxes E^ and 
F° without oscillations are taken from Ref. [0]); da/dE^ are the differential cross sections 
taken from Ref. |^, k^ is the normalized distribution of neutrino production points, h is the 



height of production, e{Ee) is the detection efficiency of the electron, "^ is the "dispersion" 
function which describes deviation of lepton zenith angle from the neutrino zenith angle ( For 
details see Ref. g]). 

Notice that the integration over the neutrino zenith angle and neutrino energy leads to 
a significant smearing of the Qi, dependence. The average angle between the neutrino and 
the outgoing charged lepton is almost 60° in the sub-GeV range and it is about 15° — 20° in 
the multi-GeV region. Neutrinos and antineutrinos of a given fiavor are not distinguished 
in the atmospheric neutrino experiments, so that the neutrino and antineutrino signals are 
summed in Eq. (|T^) which also leads to weakening of the oscillation effect. 

According to Eq. (|1^) and Eq. (|T8|) the relative excess of the e-like events, eg, can be 
represented as 

e, = ^-l^P(0^)(r(0e)4-l), (19) 

e 

where -P(Oe) is the probability averaged over appropriate energy and zenith angle intervals 
as well as over neutrinos and antineutrinos; r(0e) is the effective ratio of the electron and 
muon neutrino fluxes for a given energy and angle intervals. 

The up-down asymmetry is given by 
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Table 1: The excess integrated over Ge, eg, for different choices of Arn^i and sin^ 29 



12- 



where 



n: 



u 
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and A''e(0e) are given in Eq. (plSf). This definition includes both asymmetry of the original 
neutrino flux (in particular, due to geomagnetic effect) and the asymmetry due to oscilla- 
tions. The asymmetry due to oscillations only can be estimated as 



j^U/D 



^up ^down 
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where e^ and ef are values of the excess integrated over same zenith angle bins as in Eq. (21). 

Let us consider the sub-GeV events which correspond to the limit of integration in 
Eq. ([18|) , p < 1.33 GeV. In Fig. 1 we show the zenith angle dependences of the excess of the 
e-like events for sin^20i2 = 0.9, sin^ ^23 = 0.8 and different values of Am^i- The following 
remarks are in order. 

(i). The excess increases rapidly with Am\^ in correspondence with behaviour of the 
probability. In the first vertical upward bin the excess can reach 12.3 % for Am^i = 2 • 10-^ 
eV^. The excess integrated over all bins for different values of Am\^ is shown in the Table 1. 

(ii). The excess has rather weak zenith angle dependence. For instance, the up- 
down asymmetry equals A^l^ = 1.5 — 4.0 % for Amli > 10-^ eV^, sin^ 26'i2 = 0.9 and 



sin^^sa = 0.8. 

(ii). The excess depends very weakly on the solar mixing angle sin^ 2^12 (see Fig. 2 
and Table 1). This is related to weak dependence of sin^ 2^^ on sin^ 26'i2 in the region of 
parameters under consideration. The decrease of sin^ 26'i2 is compensated by the shift of 
the resonance to larger energies. 

(iv). The excess increases with the decrease of sin^2023- According to Eq. (|19[) the 
dependence of the excess on the mixing angle of leading channel, ^23 is determined by the 
factor (^(0e)c23 — !)• Therefore for arbitrary value of 6'23 the excess can be found using 
ee(0.8) (the calculated excess for sin^ 26'23 = 0.8; see Table 1): 

o „ , , , CnorfOp) — 1 , , cos 2^23 

e, sin^ 2^23 = e^ 0.8 ■ ^^ \'\ ^ ee 0.8 -. 23 

^ ^' ^ ' 0.723r(ee) - 1 0.45 ^ ' 

For sin^ 26'23 = 0.7, 0.9, 0.95 and r = 2 the excess (in the units of ee(0.8)) equals 1.23, 

0.716, 0.49 correspondingly. The excess disappears when the mixing in the leading channel 

approaches the maximal one. Notice, however that even for sin^ 26'23 = 0.95 an appreciable 

excess still survives. 

The excess decreases with increase of energy of the selected events. This corresponds 
to decrease of the mixing parameter in medium: sin^ 26^ oc 1/-E^ in the matter dominated 
range. In Fig. 3 and Fig. 4 we show the excess of the e-like events (as the function of the 
zenith angle Gg) for the low energy part of the sub-GeV sample with the electron momentum 
p < 0.4 GeV, and for the high energy part of the sample with p > 0.4 GeV. 

For the low energy sample (Fig. 3) the excess can reach 12 % in the upward vertical 
bin for Amli = 2 • 10"^ eV^. The excess has weak zenith angle dependence, e.g. for 
Am2i = 10~^ eV^ it decreases from 9% to 7.6 % with increases of the zenith angle . 

The integrated excess ee equals 11.3 % and 8.4 % for the two indicated values of Am|j^. 

In the high energy sample (Fig. 4) the excess in the vertical bin also reaches 12 % 
(Amji = 2 10~^ eV^), but it decreases rapidly with Arri^i. In this sample the zenith angle 
dependence is very strong, e.g. for Am2i = 10~^ eV^ the excess decreases from 5% in the 

10 



vertical upward bin to 0.7% in the vertical down going bin, the asymmetry J^[^ = 2 — 5% 
for Am2i ~ 10~"^ eV^. The integrated excess is smaller than in the low energy sample: eg = 
6% and 2.5 % for Am^^ = 2 ■ 10"^ eV^ and 10"^ eV^ respectively. 

Thus, with increase of energy the up-down asymmetry becomes more profound. However 
the absolute value of the excess decreases. 

Let us now consider the multi-GeV events. The zenith angle distributions of the excess 
for different values of Am2i are shown in Fig. 4. Let us compare the effects in the multi-GeV 
and sub-GeV samples: 

(i) Typical energy of neutrinos which produce multi - GeV events, Emc ~ 3 — 4 GeV, 
is 4 - 5 times larger than the energy, Egc, in the sub-GeV sample. Therefore the mixing 
parameter and the probability are suppressed by factor {Emcl EsgY ~ 16 — 25. 

(ii) The ratio of the fluxes, r, for the multi-GeV range is about r ^ 3 (for upward bin), 
so that the value of factor r(Oe)c23 — 1 turns out to be 2 - 3 times larger than in the sub-GeV 
sample. 

(iii) Since the average angle between the neutrino and produced charge lepton is smaller 
in the multi-GeV range, the averaging effect is smaller. In particular, in the vertical bin 
the relative contribution of trajectories which cross the core of the Earth is larger. The 
transition probability for the core crossing trajectories is slightly enhanced. 

(iv) As the result of the interplay of these factors, the excess in the multi-GeV range is 
5-7 times smaller than in the sub-GeV region. 

The zenith angle dependence of the excess is stronger. The excess decreases with | cos 0| 
and it disappears for the horizontal bin. The up-down asymmetry due to oscillations can 
reach A^l^ = 2 % for ^m\^ = 2 ■ 10~^ eV^ and it decreases rapidly with Am^i- For 
Am2i = 10~^ eV^ the excess and the asymmetry are below 1 %. 
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5 /i -like events and the Double Ratio 

In contrast to electron neutrinos, the muon neutrinos have both small and large Am'^ modes 
of oscillations. According to Eq. ([I^ ) the flux can be written as 



-^ = 1 - sin^ 2^23 sin^ ^ - e^ - e,nt, (24) 

where the first two terms correspond to the standard z/^ <-> u^ survival probability with 
slightly modified phase (see below); 



r-2 / X r^ 



e^ 



^ (rcl, - l) P2 - ^6e (25) 

rp \ / tp 



J- \ / y 



is the correction due to z/^ — z/g oscillations and 



eint = (1 - \/l - P2) sin' 2^23 cos ^ -P2 sin' 2^23 cos (f) (26) 

is the term which describes the interference of the oscillations with large and small Am? . 
The phase can be estimated in the limit of matter dominance as 

0^(Am^, + 4Am^i) — . (27) 

The first term here is the standard vacuum oscillation phase, whereas the second term is 
the correction due to 3z/ mixing. For the sub-GeV sample this correction can be neglected. 
Indeed, in the upward bins the oscillations due to Am\^ are averaged and an additional 
small contribution from the second term (associated to Ar/i'i) play no role. For downward 
bins with the average distance L ~ 100 km the second term is negligible. 

Due to strong averaging effect the correction emt can be neglected. As the result, the 
number of /i-like events can be written as 

N, = Nl\P{v,^v,)~-e,], (28) 

where P is the averaged (over energy and zenith angle) two neutrino probability. For 
maximal mixing the correction e^ is very small. With decrease of sin' 2^23 the survival 
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probability increases AP ~ l/2Asin^2^23- At the same time, the correction e^ increases 
too, thus partly compensating of AP. The compensation effect depends strongly on Am\^ 
and sin^26'23. For Aml^ = 10""^ eV^ and sin^ 26123 = 0.8 we find e^ = 2 %, whereas 
AP = 10 %. For sin^ 26^23 = 0.9 the corresponding numbers are e^ = 1.5 % and AP = 5 %. 

The double ratio R^/e = {N^/Ne){N^/N^) can be written as 

o pmax 1- 0.5 Sin^ 2^23 -e^ .^Qx 

where P^°f is the double ratio of two neutrino oscillations with maximal mixing. In the 
double ratio both corrections e^ and eg compensate the decrease sin^ 2^23- For sin^ 2^23 = 0.9 
we find AP = 5 %, e^ = 1.5 % and Cg = 3 %, so that the total increase of the double ration 
is very small: 0.5%. 

The corrections e^ and Cg are substantially smaller in the multi-GeV range, where one 

jmG ^ OsG 



would expect K^P > Rf^ 



6 Discussion and Conclusions 

We have considered the oscillation effects in the atmospheric neutrinos induced by Am^i 
and sin^ 2^2i from the region of large mixing MSW solution of the solar neutrino problem. 

The oscillations can lead to the observable excess of the e-like events in the sub-GeV 
sample with the following properties. The maximal excess is in the low energy part of the 
sub-GeV sample: the integrated excess can range from 12 % to 3% for Am^ decreasing from 
2 ■ 10^"^ to 4 ■ 10"^ eV^. The effect decreases with increase of energy of the sample. For high 
energy part of the sub-GeV range the integrated excess can reach 6 %. For the multi-GeV 
sample the excess is below 1.5 %. 

Notice that the excess of the e-like events can also be due to oscillations induced by 
the large Am'^i responsible for the leading channel of the atmospheric neutrino oscillations, 
provided that there is some admixture of the z/g - flavor in the u^ state. The two effects differ 
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by the energy dependence. The excess due to large Am^ oscillations increases with energy: 
the excess should be substantially stronger in the multi-GeV sample |Ty]. The latest SK 



data show stronger effect in the sub-GeV sample (especially in the upward going bins) thus 
preferring the solar Am'^ effect. 

The excess due to the solar Am'^ oscillations has certain zenith angle dependences: The 
up-down asymmetry due to oscillations is very weak in the low energy part of the sub-GeV 
sample, and it is strong in the high energy part of the sample. In the multi-GeV range the 
asymmetry is also strong, however the excess itself is much weaker. 

These properties will allow one to distinguish the oscillation effect from the normaliza- 
tion of fluxes (especially in future high statistics and high precision experiments). 

The excess depends on the mixing angle sin^ 2^23 responsible for the leading channel of 
oscillations. The effect decreases with sin^ 2^23 and it is strongly suppressed for maximal 
mixing, sin^ 2^23 = 1- Therefore, it is impossible to exclude the large mixing MSW solution 
from the atmospheric neutrino data or even put some bounds on parameters, unless the 
sin^ 2^23 will be measured with a good precision. On the contrary, if the excess with described 
properties will be established, this will both confirm the large mixing MSW solution and 
show that sin^ 2^23 differs from 1 . 

In this paper we discussed the excess of the e-like events. For ^23 > 45° one expects the 
suppression of the e-like events which, in fact, disfavored by the present data. Therefore, if 
Am2i will be further restricted by the solar neutrino observations, one will be able to put 
the bound on sin^ 26*23. 

Comparison of the predicted excess with data (Fig. 1-5) shows reasonable agreement. 
For Am2i ~ 10^'* eV^ the excess due to oscillations reproduces both the size and the zenith 
angle distribution of the observed excess in the sub-Gev range. It also gives reasonable 
description of the data in the low energy part of the sub-GeV sample. At the same time, 
the predicted excess is smaller than the detected one in the high energy sub-GeV range, 
and also in the multi-GeV range. Notice, however, that in the high energy sub-GeV range, 
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and especially in multi-GeV range, the effects of the large A?ti|^ (neglected here) can be 
important [|T^ . They can change the signal by 15 % in the multi-GeV sample and by 5 - 6 % 
in the high energy part of the sub-GeV sample. Also non-trivial interference of small and 
large Am^ oscillation effects is possible pO[. It is clear however, that with present statistics 
it is impossible to make definite conclusions. In fact, it could be that the explanation of 
the data will require some interplay of the normalization of fluxes and different oscillation 
effects. 

The excess has the positive up-down asymmetry. The positive asymmetry at low en- 
ergies can be reconciled with the negative asymmetry at high energies (multi-GeV sample) 
indicated by the present data. Indeed, the excess due to solar Am^ (with positive asymme- 
try) decreases with increase of energy and in the multi-GeV range the dominant effect will 
be due to oscillations with large Am^ and nonzero ^13. In this case the sign of the excess is 



determined by factor |[T0| : 



(rsl, - 1) (30) 

(with S23 instead of C23, see Eq. (0)). In the multi-GeV range the average ratio of the fluxes 
r ^ 2.5, and for S23 = 0.277 we get negative value of the factor (|30D : —0.3. Indeed, the latest 



SK data indicate a negative up-down asymmetry at high energies. The interpretation of the 
multi-GeV zenith angle distribution can require some interplay of the normalization which 
explains the excess of events in the down-going bins (cosO = 0.2 -^ 1) and the oscillations 
which suppress the number of events in the upward going bins (cos© = —1 -. — 0.2). 

Further studies of solar neutrinos (in particular, searches for recoil electron spectrum 
distortion and the day-night effect) will allow one to prove or disprove the large mixing 
MSW solution. Observations of signiflcant day-night effect will imply small values of Arn^^ 
and therefore small excess of the e-like events. On the contrary, weak day-night effect and 
upturn of the distortion of the recoil electron spectrum at low energies will testify for large 
Am|^ and large excess. The KAMLAND experiment [^ will test whole the range of the 
large mixing angle solution. 
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If the effects discussed in this paper will be confirmed we can be left with the "bi-large 
mixing scheme" with large (but not maximal) mixing between neighboring generations and 
small e — r mixing. The mass hierarchy will be rather weak: m2/m^ ~ 0.1 — 0.3 which will 
allow one to explain large mixing without special arrangements. 

Note added 

After this work had been accomplished the paper pBf has appeared in which the im- 



plications of the low-energy SK data for LMA solution of the solar neutrino problem are 
discussed. It is claimed that oscillations lead always to decrease of number of e-like events 
in comparison with no oscillation case. This statement contradicts our results. The error 
follows from using the ratio of numbers of events, R^S-, instead of ratio of the neutrino 



fluxes (r in our notation) in formulas (2.22)-(3.9) of Ref. [Q. Also matter effects must be 
included. 
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Figure captions 

Fig. 1. Zenith angle distribution of the excess of the e-hke events in the sub-GeV range for 
sin^ 26*23 = 0.8, sin^ 26*12 = 0.9 and different values Amgi- The points are the 735 days data 
of the Super-Kamiokande @]. 

Fig. 2. Zenith angle distribution of the excess of the e-like events in the sub-GeV range 
for different values of sin^26^i2 and for two values of Am2i = 4 10~^ eV^ and 2 10~^ eV^; 
sin^ 26*23 = 0.8. The points are the 735 days data of the Super-Kamiokande 0. 

Fig. 3. The same as in Fig.l for the low energy part of the sub-GeV range {p < 0.4 GeV). 

Fig. 4. The same as in Fig.l for the high energy part of the sub-GeV range {p > 0.4 GeV). 

Fig. 5. The same as in Fig.l for the multi-GeV range. 
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